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I. Abstract 

Results of five years of Earth-orbital exposure on mechanical properties of glasses 
indicate that radiation effects on mechanical properties of glasses, for the glasses 
examined, are less than the probable error of measurement. During the 5 year exposure, 
7 micrometeorite or space debris impacts occurred on the samples examined. These im- 
pacts were located in locations which were not subjected to effective mechanical testing 
hence limited information on their influence upon mechanical strength was obtained. 
Combination of these results with micrometeorite and space debris impact frequency 
obtained by other experiments permits estimates of the failure probability of glasses 
exposed to mechanical loading under earth-orbit conditions. This probabilistic failure 
prediction is described and illustrated with examples. 

II. Introduction 

Design of glass components employed in earth orbit has a long and rich history 1 and 
glass elements of orbital systems have not contributed to any known system failures. 
Growth of orbital debris and reports of micrometeorite impacts on space systems have 
heightened concern for safe design of glass components of orbital systems. To the 
authors knowledge no analytical techniques are currently available to predict failure 
probabilities of glass components in orbital systems. 

Extensive data is available from the NASA LDEF (Long Duration Exposure Facility) 
and other flights to permit characterization of the size and frequency of micrometeorite 
and space debris impacts on orbital systems as is shown in Figure 1. The authors 
recently completed 2 analyses of mechanical strengths of a group of glass samples 
exposed on LDEF. The total number of impact events suffered by these glasses was 
small and the data obtained are insufficient to permit statistically valid predictions of 
failure probabilities of glass in space. This paper combines the results of careful analysis 
of a small group of impact events on glasses with the larger collection of micrometeorite 
and debris impact data to produce a prediction of the probability of failure of a glass 
exposed in earth orbit as a function of stress on the glass. 
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Ideally the determination of failure probability of a stressed glass sample in earth orbit 
would be determined by exposure of a large group of stressed samples for a period of 
time sufficient to obtain statistically valid exposures to the impact of a large number of 
micrometeorites or space debris. Unfortunately such experiments would require a large 
number of mechanical loading fixtures and associated instrumentation to ascertain the 
failure times for a variety of stress levels. Such an experiment is beyond the capability of 
existing hardware. 

The experiment which has been completed is the exposure of 120 glass and glass ceramic 
samples in a passive unloaded state for a period of 5.8 years in earth orbit as previously 
described. These experiments resulted in the conclusion that the mechanical strength 
and optical transmission of those samples was, within the limitations of the measurement 
and experiment, unchanged by earth-orbit exposure. The samples did however suffer 
seven micrometeorite or debris impacts which clearly degraded the mechanical strength 
at the impact location. These impacts, because of their statistically determined locations, 
did not occur at locations permitting proper determination of their effect on the glass 
mechanical properties. 

It was possible to determine the maximum stress concentration for several impacts, and 
this data provides the basis of the calculations reported here. The conclusion of this 
analysis 3 was that the damage field associated with the impacts had a depth to diameter 
of crater of less than 1/4 which resulted in a stress concentration of not more than 2. 
Data gathered by other LDEF investigators 4 on the flux of particles sharing the orbit 
with spacecraft was incorporated into the mechanical damage calculations. The results 
enable prediction of the life expectancy of glass components in space. The design of 
glass components for earth orbiters may be refined by application of this result. 


m. Experimental Procedure 

Samples, from which the strength measurement results 2 (Table I) of glasses 
exposed on the LDEF (Table II) were obtained, were examined for impact sites from 
micrometeorite or debris. The stress field during biaxial flexure testing 5 was related to 
the position of the micrometeorite or debris impacts. An estimate of the maximum 
stress at the site of the impact was calculated from the failure strength of the sample 
containing the impact. Fracture mechanics 6 were applied using the estimated stresses to 
determine the maximum flaw size capable of withstanding the stresses before failure. 

The crack depths determined from the relation between critical flaw size and applied 
stress, were compared to radial crack lengths and crater diameters measured from 
scanning electron micrographs of the impacts. Particle size was calculated from crater 
size using cratering mechanics 7 *. Statistical information on particle size and flux 4 was 
incorporated into the probability of failure calculation. A model of a circular window 
exposed to conditions in a low earth-orbit environment was developed. The stresses on 
the surface of an 80 cm diameter window under a pressure of 1 atmosphere were deter- 
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mined. Equations were developed relating the stress surrounding an impact to the flux 
of particles in a low earth-orbit. These equations were integrated with respect to the 
dimensions of the window to determine the expected time to failure. Probability of 
failure was calculated using the Poisson distribution. 

IV. Results 

Stresses at the point of the micrometeorite or space debris impacts were found for each 
of the impacted samples described in Table I. Table III reports observed crater 
diameters, spall diameters and melt pit diameters. Particle size (D p ) was computed from 
crater diameter (dj with the aid of defined cratering relationships. Cratering mechanics 
provides established relationships between projectile size and crater size for projectiles 
with similar densities and velocities. A relation was found between radial cracking, melt 
pit diameter, and crater diameter for the impacted LDEF glasses. According to the data 
in Table III, radial crack lengths are on the order of five times the melt pit diameter, 
and between two and three times the crater diameter. Through this series of relation- 
ships the depth of the surface flaw created during impact in glass may be estimated for 
any given particle size assuming the particles have similar densities and velocities. 

Particle flux versus particle size data reported by Kinard 4 as illustrated in Figure 1 was 
integrated with the cratering data for the 7 impacts. The surface damage per unit time 
per unit area leading to failure was determined for a glass component in the earth-orbit 
environment. The expected time to failure depends on the surface stresses of the 
component and the depth of impact damage. 

Failure of glass occurs when the stress surrounding a flaw exceeds the strength of the 
glass. Maximum tensile stresses due to bending occur at the surface where flaws 
intensify local stresses. The flaw depth and crack tip radius determine the stress 
concentration associated with the flaw. A flaw created by the impact of a 
micrometeorite or space debris particle with an orbiter window will cause failure if the 
stress concentration factor at the impact site and the stress applied at that location 
exceeds the strength of the glass. Particles producing flaws with depths less than the 
critical crack length will not cause failure, thus failure from impact depends on flaw 
depth and local stress. Local surface stress is determined by the applied pressure and 
window geometry. Particles producing flaws larger than the critical size will cause failure 
in glass, thus the threshold particle size depends on the specific stress of each surface 
element. The total number of impacts leading to failure for a window depends on the 
product of the area under constant stress, the time of exposure, and the flux of particles 
exceeding the threshold size to create critical flaws. The average cumulative flux ( <p ) of 
micrometeorites and space debris arriving over a period of one year, per square meter of 
area was found from Kinard’s data to fit a line approximated by: 

<(> = bD a ( 1 ) 

where b = 4 x 10 - ®, m = 2.1, and D is the particle size. 


3 



The critical particle size (D c ) is the particle size which will produce a flaw with a 
depth exceeding the critical crack length. D in equation (1) may be replaced by D c to 
yield an equation which describes the flux of particles causing failure. The critical 
particle size was determined from cratering relations and the critical flaw size. The 
critical flaw size (a e ) was found by applying the stress concentration equation for a 
semielliptical surface flaw which is given as: 



( 2 ) 


where is the theoretical strength, a. is the applied stress, a is the crack depth, and p is 
the curvature of the crack tip. Failure occurs when the crack depth equals the critical 
flaw size needed to increase the stress at the crack tip to the theoretical strength of the 
glass. This occurs in glasses impacted by micrometeorites or space debris when the 
particle size and velocity is great enough to produce damage greater than or equal to the 
critical flaw size. The average velocity of particles for the trailing-side of LDEF was 
estimated 9 to be on the order of 13 km/s. Flaw depth may be determined from projec- 
tile size, assuming the velocity of impacting particles for the samples in this investigation 
is constant. The critical particle diameter (D c ) may be substituted in the stress concentra- 
tion equation along with the radius of curvature for the median cracks. The radius of 
curvature was calculated by the authors 3 to be about 4 x 10' 3 yxm. The equation as 
applied to the impacted glasses becomes: 


°c 





(3) 


which yields the critical stress leading to failure for a given particle size. The particle 
size creating a critical flaw is therefore: 


D c 


= ( 0ai ) 2 SL 
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Equation (1), written in terms of the critical particle size, defines the cumulative particle 
flux causing failure in earth-orbit glasses. Substituting equation (4) for D c in equation 
(1) yields a flux equation, which is dependent on the applied stress, given by: 


♦ 


(o) 
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The expected number of impacts causing failure per unit time per unit area was found by 
applying the stress-dependent flux equation over the surface of the glass. The calcula- 
tion is the product of <p w and the area under constant stress. Generally, glasses used in 
space applications would not be expected to have constant stress over the entire surface 
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exposed to the space environment; hence, the total number of expected fatal impacts per 
unit time is found by integrating the flux equation for particles greater than or equal to 
the critical size over each area increment under critical stress. The expected number of 
impacts leading to failure was determined by integrating over the surface: 

jj<|> (0) dadA (6) 


Application to an Orbiter Window 


The stress distribution on the surface of a glass window with circular geometry is deter- 
mined by 10 : 



6 P (3+v) 
16 t 2 


(R 2 -z 2 ) 


(7) 


where P is the applied pressure, R is the outer radius, r is the radius, and v is Poisson’s 
ratio, which for silicate glasses is about 0.3. The maximum stress occurs in the geometric 
center of a simply supported circular window. The stress decreases as r approaches the 
outer radius and is constant in any annulus. The area of each annulus with constant 
width increases as r increases; thus area limits the number of particles causing failure in 
the near geometric center of the window, and stress, which tends to zero, limits the 
number of particles causing failure in the outer annuli. This premise is valid for particle 
sizes less than those required for penetration, which limits the flux of critical particles to 
be greater than the penetration flux. The flux of critical impacts is illustrated by a graph 
of the number of impacts causing failure per annulus per year as a function of radius. 
This function is represented by Figure 2 for a window under 1 atmosphere of pressure 
with a diameter of 0.8 meters (area = 0,5 m 2 ) and a thickness of 3 cm. The area under 
the curve is the total number of impacts which the window is expected to receive per 
year. The surface integral giving the total number of expected fatal impacts for the 
circular window is: 


2nr x 

f zdrdd = 2rcj*$ (r) rdx (8) 

0 0 0 

where $ (r) is a function of the cumulative particle flux , surface tensile stress, and annulus 
area, the latter two being dependent on the radius r. The value of this integral is well 
approximated by taking the sum of the impacts causing failure per year per annulus over 
all annuli by: 
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The number of years expected before failure is the inverse of this sum. 


The probability of failure may be determined from the cumulative number of impacts 
causing failure. The impact events may be considered as Poisson arrivals over a given 
period of time. The Poisson probability function is given by 11 : 


P (x,l) 


g~U x 
x l 


( 10 ) 


where X is the sum found by equation (9) and x is the Poisson variable equal to 0 if no 
failure occurs. The probability that failure will occur in time t is given by: 


P M«u„ ■ 1 ' 


( 11 ) 


Figure 3 shows the expected years to failure as a function of applied pressure for several 
different thicknesses. Windows with no applied stress are limited by the flux of particles 
large enough to cause penetration. The expected time for one impact causing failure in 
glass with a thickness of 3 cm, a diameter of 0.8 meters, under a pressure of 1 atmo- 
sphere is 35 years, which is a reasonable requirement for earth-orbit space station 
components. Equation (11) estimates a 0.84 probability that the glass would fail in 35 
years and a 5% chance of failure in one year. A window with similar geometry requires 
a thickness over 6.5 cm to achieve a 0.008 failure probability in 35 years. A 4.7 cm thick 
window with no applied stress yields an equivalent failure probability. 


V. Discussion 

Reliance upon a single LAYER window does not permit satisfaction of design lifetime 
criteria for earth orbiters. The probability 12 of success for components affecting crew 
safety in the current shuttle program is currently 0.992. The failure probability of a 
window may be significantly decreased by using multi-layer windows. In fact, current 
window design for spacecraft in the shuttle program employs multi-layer fused-silica and 
alumino-silicate windows 13 . Failure of the window is dependent on the joint probability 
of each component. An independent joint probability assumes that once the integrity of 
an outer layer is violated then the outer window offers no protection to inner windows. 
This assumption produces a conservative estimate of the failure probability for multi- 
layer windows. The outer glass in a multi-layer window requires a larger particle size to 
initiate failure as compared to the inner pane because the outer pane is not pressurized. 
This circumstance eliminates the dependence of failure on stress. Failure of the 
unstressed window due to micrometeorite or space debris interaction becomes dependent 
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on penetration of an impacting particle. Assuming the ratio of crack depth to particle 
size remains constant, the integrity of the window would be compromised for particles of 
1/4 the glass thickness according to cratering analysis. This assumption is valid only if 
the velocity of particles is no greater than the velocity of those particles impacting the 
trailing edge of LDEF. Cratering calculations 7 suggest penetration is dependent on 
impact energy; thus, smaller particles penetrate at higher velocities. Multi-wall protec- 
tion of critical parts in space applications is a well known means of protection and has 
been well investigated 1411 ' 16 . The outer plate is an important component since its purpose 
is to absorb the initial shock and fragment or contain the impacting particle. A failed 
outer layer not leading to immediate inner layer failure may be replaced, which reduces 
the chance that a particle will impact the inner window. 

Probability of failure for a two layer window may be determined by applying the 
probability density function of the Poisson distribution to the joint probability of failure 
of two layers. The probability is calculated by integrating the probability density 
functions over time of service. The probability density function for the windows is deter- 
mined from the integral of the joint probability function given as: 

t 

e'* 1 " \ 2 dT (12) 


where A., and X 2 are the flux of critical impacts for the inner and outer window respective- 
ly, t is the time to outer layer failure and t is the time to inner layer failure. To achieve 
a 0.992 probability of success in a two-layer window with an inner 3 cm thick window, 
the outer window would need to be 3.5 cm thick. 

The Poisson probability function is convenient for approximating binomial functions 
which depend on several variables. The probability of failure found from the Poisson 
distribution approximates the value determined directly from the flux equation of the 
expected number of impacts leading to failure for less than 0.1 expected fatal impacts 
per year. For rates of more than 0.1 expected fatal impacts per year, the expected 
number of impacts is not valid for approximating the probability of a component failure. 
The Poisson limit accounts for survival chances at high probabilities of failure and is an 
accurate model when the total number of events is not defined; hence, this distribution is 
widely used to model "real world" phenomena. 

Damage due to micrometeorite or space debris in the glasses incorporated in this study 
occurred at velocities on the order of 13 km/s. Mandeville 17 has found that craters 
formed in glass in the hyper-velocity range to be primarily dependent on velocity. 
Dependence on projectile and target density was secondary. His findings determined the 
central pit diameter was on the order of 2.5 times the particle size. Estimations of crater 
diameters determined by explosive cratering studies 18,19-20 have been determined by 
applying Pi-group scaling laws. The calculations were developed for cratering in rock 
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and soil, but were found by this author to approximate microcrater sizes on the order of 
the spall diameters listed in Table I. Assuming the scaling laws are valid for 
micrometeorite and space debris impacts over the range of particle velocities in low 
earth-orbit, flaw depths in glass may be estimated for impacts above 13 km/s using 
cratering dimensions and particle flux. Scaling laws predict that crater size is dependent 
on E 1/3 , where E is impact energy; thus, the damage from micrometeorite and space 
debris would be expected to increase as v 273 ; however, this dependence has not been con- 
firmed in earth-orbit glasses. Melting 7 occurs in glass at impact velocities greater than 
about 10 km/s. Depth of damage may be reduced by melting and viscous flow since the 
dissipated energy may relieve stresses related to crack growth. 

The velocity of particles impacting glass on an earth orbiter depends on the orientation 
of the glass with respect to the ram direction of the orbiter. The velocity 21 of particles 
impacting the ram side of earth orbiters has been estimated to be about 20 km/s; 
therefore, the flaw depths produced by micrometeorite and space debris should be a 
maximum at this orientation. The cratering process is dependent on the particle velocity; 
thus, the magnitude of melting and cratering produced by particles at velocities greater 
than 13 km/s would be expected to scale accordingly. Assuming the relation between 
crater size and depth of damage in the glasses is comparable for higher velocity impacts 
and that the flux of particles is constant, the probability of failure for glasses exposed to 
higher velocity particles may be calculated. 

Windows of an earth orbiting space-station would be susceptible to hazards of the earth 
orbit environment. Glasses integrated into such systems must be designed to maximize 
the safety of the crew members for the duration of the orbiter’s flight. Proper window 
design must include the probability of a penetrative meteorite or space debris impact 
with the window systems. This study has produced a means of predicting the lifetime 
expectancy of such windows by applying cratering relationships and fracture mechanics to 
glasses exposed in the earth orbit environment. Further studies should examine impact 
parameters determined from higher velocity impacts in order to investigate crack 
morphologies of the higher energy impacts. Information on flaws in glass produced at 
higher velocity would be useful in representing particle impacts in glass components on 
the ram side of earth orbiting modules. 

VI. Conclusions 

Measured micrometeorite and debris flux, and particle size may be used to estimate 
failure probabilities of glasses in earth-orbit. 

Failure probability of windows depends on design, pressure and geometry. The use of a 
multi-pane design significantly reduces the probability of window failure. 

Earth orbit life expectancy of glasses may be computed from particle flux and cratering 
relationships. 
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Table I. Location of Impacts with respect to Stress Contours 


Sample 

Stress 

Type 

Stress 

Contour 

Failure Stress 
(MPa) 


Minimum 

(MPa) 

Fused Silica 

0.4 

103 

39 


Zerodur 

0.2 

125 

31 


Zerodur 

0.5 

125 

62 


Zerodur 

0.3 

135 

41 


Pyrex 

0.1 

106 

11 


Soda-Lime 

0.4 

105 

37 


BK-7 


130 

13 
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Table in. Crater Dimensions for Glasses Impacted by Micrometeorites or Space 
Debris 



Melt 

Crater 

Sp; 


Diameter 

Diameter 

Diameter 


(/xm) 

(urn) 

(Mm) 

BK-7 

40 

100 

200 

Fused Silica 

50 

120 

250 

Soda-Lime-Silica 

80 

175 

475 

Pyrex 

85 

200 

400 

Zerodur* 

75 

200 

400 

Zerodur* 

50 

150 

300 

Zerodur** 

No Melt 

100 

275 


* This sample received two impacts. 
** This sample exhibited no melting. 
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